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ABSTRACT 



How should you choose a good set of (say) 48 planes in four dimensions? More generally, 

how do you find packings in Grassmannian spaces? In this article I give a brief introduction 

to the work that I have been doing on this problem in collaboration with A. R. Calderbank, J 

H. Conway, R. H. Hardin, E. M. Rains and P. W. Shor. We have found many nice examples of 

specific packings (70 4-spaces in 8-space, for instance), several general constructions, and an 

| embedding theorem which shows that a packing in Grassmannian space G(m, n) is a subset of 

■ a sphere in R , D = (m + 2)(m — l)/2, and leads to a proof that many of our packings are 

00 , 

' optimal. There are a number of interesting unsolved problems. 

(N ' 
O 

_> I 1. Introduction 

In my talk at the Yamagata conference on "Algebraic Combinatorics and Related Topics" 
(November 1997) I discussed two problems, (a) finding packings in Grassmannian manifolds, 
^ and (b) constructing error-correcting codes for quantum computation, and tried to show how 

our purely numerical investigations into the first problem had led to theoretical advances in 



both subjects. This work has been presented in a series of papers, [10], pJ] , ]i~2"[ , |i~3| ], [37]. In 
the present paper I will give a brief introduction to our work on the first problem, referring 
the reader to the above references for further information. 

For a long time I have been interested in various kinds of packing problems: packings in 
Hamming space (i.e. error-correcting codes p0|), in Euclidean space (i.e. the sphere-packing 



problem |l4|| ), or on the sphere (i.e. spherical codes [14, Chap. 3], |23f| ). R. H. Hardin 
and I have had some success in using numerical optimization techniques to search for spherical 
codes and other kinds of geometrical designs. This work is described in the papers [20|, [21|, 



p4}| , |^5| , |2(| , [^?J (others are in preparation) . Figure 1 (on the next page) shows an example 
that arose from our search for good coverings, i.e. sets of points with small covering radius. 
The figure shows 48002 points on the sphere with average angular separation of about 1 degree 
(the minimal angular separation is .802 degrees, the maximal angle is 1.077, the average angle 
is .99948 and the standard deviation is .047). 

A few years ago, a statistician, Dianne Cook, asked if we could apply the same techniques 
to find packings of (say) 48 planes in 4-dimensional Euclidean space Mr. More generally, for 
given values of N, n and m, how should one arrange N n-dimensional subspaces of M m (all 
passing through the origin) so that they are as far apart as possible? 

The rest of the paper will describe what happened when we tried to attack this problem. 

2. Packings in Grassmannian Space 

The Grassmannian space G(rn, n) is the set of all n-dimensional subspaces of real Euclidean 
m-dimensional space M m . This is a homogeneous space isomorphic to 0{m) / (0(n) x 0(m—n)), 
and forms a compact Riemannian manifold of dimension n{m — n). 

Before we can investigate packings in this space, we must decide how to measure the 
distance between two n-spaces P, Q £ G(m,n). The principal angles 9i,...,9 n £ [0,7r/2] 
between P and Q are defined by (we follow Jl9| , p. 584) 

cos 0i = max max u ■ v = u% • Uj , 
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for i = 1, . . . , n, subject to u ■ u = v ■ v = 1, u ■ Uj = 0, v ■ Vj = (1 < j < i — 1). The vectors 
{iii} and {vj} are principal vectors corresponding to the pair P and Q. 

Wong |42j shows that the geodesic distance on G(m, n) between P and Q isQ 

d g (P,Q) = j0l + .~ + 0l. (1) 

However, this definition has one drawback: it is not everywhere differ entiable. Consider the 
case n = 1, for example, and hold one line P fixed while rotating another line Q (both passing 
through the origin). As the angle <j> between P and Q increases from to ir, the principal angle 
6\ increases from to tt/2 and then falls to 0, and is non-differentiable at tt/2 (see Fig. ||). 




tt/2 (f) 



Figure 2: Principal angle 0\ between two lines as the angle between them increases from to 

7T. 

Although one might expect this non-differentiability to be a mere technicality, it does in 
fact cause considerable difficulties for our optimizer, especially in higher dimensions in cases 
when many distances fall in the neighborhood of singular points of d g . 

An alternative measure of distance, which we call the chordal distance, is given by 



d c {P,Q) = Vsin i 0i + ••• + sin 2 ^ n . (2) 

The reason for the name will be revealed later. This approximates the geodesic distance when 
the subspaces are close, has the property that its square is differentiable everywhere, and, as 
we shall attempt to demonstrate, has a number of other desirable features. 

A third definition has been used by Asimov P] and Golub and Van Loan [19], p. 584, 
namely 

d m (P,Q) = max 0i . 

i=l,...,n 

This shares the vices of the geodesic distance. 

1 The geodesic distance is unique except for the single case of G(4, 2). 
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Of course for n = 1 all three definitions are equivalent, in the sense that they lead to the 
same optimal packings. 

We can now state the packing problem: given N, n, m, find a set of re-spaces Pi, ... , P/v £ 
G(m,n) so that min d(Pi,Pj) is as large as possible, where d is either geodesic or chordal dis- 
tance. Since G(m, n) is compact, the problem is well-defined. Because G(m, n) and G(m, m—n) 
are essentially the same space, we may assume n < to/2. 

We also need some further terminology. A generator matrix for an n-space P £ G(m, n) is 
an n x m matrix whose rows span P. The orthogonal group 0(m) acts on G(m,n) by right 
multiplication of generator matrices. The automorphism group of a subset {P\, . . . ,Pn} C 
G(m,n) is the subset of 0(m) which fixes or permutes these planes. 

By applying a suitable element of 0(m) and choosing appropriate basis vectors for the 
spaces, any given pair of n-spaces P, Q with n < to/2 can be assumed to have generator 
matrices 

"1 ••• ••• ••• o- 

(3) 
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respectively, where 9\, . . . , 9 n are the principal angles Theorem 2) 

3. Packing lines in 



The next step was to apply a suitably modified version of our optimizer in order to find 
a collection of good packings. The basic algorithm is described in and we shall not say 
much about it here. We had to overcome various obstacles caused by the nondifferentiability 
of the geodesic distance. Whereas we could compute partial derivatives of the chordal distance 
analytically, for the geodesic distance we used numerical differentiation. 

Our initial computations concerned packings in G(m, 1); that is, packings of lines through 
the origin in R m , or equivalently the problem of packing points on the sphere in R m so that 
if P is a point in the packing, so is —P. Such sets of points form antipodal spherical codes. 
This is a classical problem, of course, that is the subject of many papers. The 3-dimensional 



problem is discussed for example in Qj, f32|, g§, @, g§, @], |3g], @0|. 
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Nevertheless we were able to find many packings that were better than those previously 



known, even in three dimensions. We refer the reader to [13], [22] for the details. Here we 
will just give one table, showing the minimal angle 6\ of the best packing we have found of 
N lines, and for comparison the minimal angle 6'i of the best packing known of 2N points on 
S 2 (taken from p2[ ). We see that requiring a packing of 2N points on S 2 to be antipodal is 
a definite handicap: only in the cases of 6 and 12 points do the antipodal and unrestricted 
packings coincide. Decimals in the tables have been rounded to four decimal places. 

The last two columns of the table specify the largest automorphism groupf] we have found 
of any such best antipodal packing of 2N points. The fourth column gives the order of the 
group and its name both in the orbifold notation (cf. |i~5]| ) and as the double cover of a rotation 
group. The symbol ±G indicates that the group consists of the matrices ±M for M £ G, where 
G is a cyclic (C), dihedral (£>), tetrahedral (T), octahedral (O) or icosahedral (I) group. In 
each case the subscript gives the order of the rotation group. In some cases the best packings 
can be obtained by taking the diameters of a known polyhedron, and if so this is indicated in 
the final column of the table. 

The entries for N < 6 were shown to be optimal by Fejes Toth in 1965 |l^] (see also 
Rosenfeld |33|), and the 7-line arrangement will be proved optimal in Section 5. The solutions 
for N > 8 are the best found with over 15000 random starts with our optimizer. There is no 
guarantee that these are optimal, but experience with similar problems suggests that they will 
be hard to beat and in any case will be not far from optimal. 

For N = 1,2,3,6 the solutions are known to be unique, for N = 4 there are precisely two 
solutions (|]18|, |3|], H), and for N = 5,7,8 the solutions appear to be unique. For larger 
values of N, however, the solutions are often not unique. For N = 9 there are two different 
solutions, and in the range N < 30 the solutions for 10, 22, 25, 27, 29 lines (and possibly 
others) contain lines that "rattle" , that is, lines which can be moved freely over a small range 
of angles without affecting the minimal angle. The table only goes as far as N = 28 lines. For 



more a extensive table in G(3, 1) and for tables of packings in G(m, 1) for m > 3, see [13]. 

4. Packing planes in r 4 

We were naturally very interested to see what would happen when we studied packings 

in Cr(4, 2), that is, packings of planes in Mr. Here for the first time we felt we were sailing 
2 That is, the subgroup of 0(3) that fixes or permutes the 2N points. 
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Table 1: Best packings known of N lines in R 3 {9\ gives angle). 



N 


min 9\ 


min 9[ 


group 


polyhedron 


2 


90.0000 


109.4712 


*224 


= ±Va 


square 


3 


90.0000 


90.0000 


*432 


= ±OlA 


octahedron 


4 


70.5288 


74.8585 


*432 


= ±OlA 


cube 


5 


63.4349 


66.1468 


2*5 


=±£>m 


pentagonal antiprism 


6 


63.4349 


63.4349 


*532 


- 1 — ^OU 


icosahedron 


7 


54.7356 


55.6706 


*432 


= ±OoA 


rhombic dodecahedron 


8 


49.6399 


52.2444 


x 






9 


47.9821 


49.5567 


3x 






9 


47.9821 


49.5567 


2* 


= ±Co 




10 


46.6746 


47.4310 


*226 


= ±£>19 


hexakis bi-antiprism 


11 


44.4031 


44.7402 


2*5 


=±Pm 




12 


41.8820 


43.6908 


*432 


= ±Oo 4 


rhombicuboctahedron 


13 


39.8131 


41.0377 


2* 






14 


38.6824 


39.3551 


x 


=±Ci 

- J -*- , l 




15 


38.1349 


38.5971 


2*5 


=±Pm 

- J - i -^ ±u 




16 


37.3774 


37.4752 


*532 


= ±J fin 


pentakis dodecahedron 


17 


35.2353 


35.8078 


x 


=±Ci 




18 


34.4088 


35.1897 


3x 






19 


33.2115 


34.2507 


X 


=±Ci 




20 


32.7071 


33.1584 


*222 


= ±V 4 




21 


32.2161 


32.5064 


5x 


=±c 5 




22 


31.8963 


31.9834 


2*3 


=±v 6 




23 


30.5062 


30.9592 


X 


=±Ci 




24 


30.1628 


30.7628 


3*2 


=±Ti2 




25 


29.2486 


29.7530 


3x 


=±c 3 




26 


28.7126 


29.1948 


2* 


=±c 2 




27 


28.2495 


28.7169 


X 


=±Ci 




28 


27.8473 


28.1480 


X 


=±Ci 





in waters where no one had been before. (Of course, as already mentioned in Section 2, we 
expected that G(4, 2) would be special.) 

Using our optimizer, we looked for the best packings we could find of N planes in R 4 (i.e. 
N points in G(4, 2)), for values of N up to about 50. 

The coordinates of the planes as found by the computer are with respect to a random 

coordinate frame, and so must be "beautified" by hand. Initially they look like this 

-0.4909573575989161 -0.5698930951299707 -0.6378400314190236 0.1653566674254520 
-0.7841182881275834 0.0542768735416466 0.4433546762436948 -0.4308702383261817 
-0.0989775328964588 -0.2213935599532070 -0.9654594563867928 -0.0952700250146647 



That is, they have no obvious structure. The reason they are called "planes" is that they are 
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indeed "plain" (at this point I held up a sheet of perfectly blank white paper; the audience 
laughed). The fun comes in trying to understand the computer output. 

For N = 2 planes, the best packings for both definitions of distance are the same: take two 
orthogonal planes, say 

+ 000 00 + 

0+00 000+' 

with principal angles w/2, so that tt/2, d\ = 2, dz = ir 2 /2. (We abbreviate +1 and — 1 by + 
and — respectively.) 

For N = 3 planes there are different answers for the two distances. For the geodesic distance 
the best packing consists (for example) of the planes 

+ -00 +0-0 +00- 

000 + ' 0+00' 00 + 0' 

and has dg = 57r 2 /18, d 2 = 1.25. For chordal distance the best packing consists (for example) 
of the planes 

10 r .5 1 -r .5 10 0-1 

1 -.5 r ' 1 -.5 -r ' 11 ' 

and has d 2 = 1.5, dg = 2tt 2 /9. A more geometrical description of these two arrangements will 
be given below. 

I postpone discussion of N = 4 and 5, and consider N = 6 next. Here the same arrangement 
appeared to be optimal for both distances. In this packing the principal angles between any 
two of the six planes are tt/2 and arccos 2/\/5, so the 6 planes lie at the vertices of a regular 
simplex in G(4, 2). This suggested that the arrangement should somehow be related to the 
icosahedron, since the angles between any two of the six diameters of the icosahedron are equal 
(i.e. these six diameters form a regular simplex in G(3, 1)). We soon realized that there is a 
classical theorem which explains this, and gives a simple way of describing G(4, 2). 

We remind the reader that any element a of SO (4) may be represented as 

a : x \- 



where x = XQ + x\i + X2j + x%k represents a point on S 3 and £, r are unit quaternions [17]. The 
pair —£, —r represent the same a. The correspondence between a and ±(l,r) is one-to-one. 

Given a plane P E £7(4, 2), let a be the element of SO (4) that fixes P and negates the points 
of the orthogonal plane P ± . Then a 2 = 1, and for this a, it is easy to see that t = ^ii+^2j+^3^ 
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and r = r\i + r^j + r^k are purely imaginary unit quaternions. This establishes the following 
result (which can be found for example in Leichtweiss |29f|). 

Theorem 1. A plane P £ G(4, 2) is represented by a pair (£,r) 6 S 2 x S 2 , with (— r) 
representing the same plane. 

Given two planes P,Q£ G(4, 2), represented by ±(£, r), ±(£' , r'), respectively, the principal 
angles 9\, #2 between them may be found as follows. Let <p (resp. ip) be the angle between £ 
and £' (resp. r and r'), with < <p, ip < tt. If <p + ip > tt, replace (p by tt — <p and ip by ir — ip, 
so that < <p + ip < tt, with (p <ip (say). Then 

h,02 = ^^, 4(P,Q) = ^^, ^(P,Q) = 1-cos^ cos0. 

A set S = {Pi, . . . ,Pn} ^ G(4, 2) is thus represented by a "binocular code" consisting of 
a set of pairs ±(£j, r^) £ 5" 2 x S 2 . We call the list of 2iV points ±£j (they need not be distinct) 
the "left code" corresponding to S, and the points ±rj the "right code". Conversely, given 
two multisets L C S 2 , R C 5 2 , each of size 2iV and closed under negation, and a bijection or 
"matching" / between them that satisfies f(—£) = —f{£), £ £ L, we obtain a set of ./V planes 
in G(4,2). 

The binocular codes for the d c -optimal packings of N = 2, . . . , 6 planes are shown in Figs. [3|, 
|i|. Except for iV = 3, the left and right codes are identical. Matching points from the left and 
right codes are labeled with the same symbol. For N < 5 the points lie in the equatorial plane, 
and for N < 4 there are repeated points. The points lie on regular figures, except for N = 4 
where the points are ±(1,0,0), ± ("^g, ^\[\-> o) • 

For N = 6 the left and right codes consist of the 12 vertices of an icosahedron (Fig. |^). Let 
these be the points 

A(0,±1,±t), A(±t,0,±1), A(±1,±t,0) , 

where A = l/\A~ + 2. The matching is obtained by mapping each point to its algebraic 
conjugate (i.e. replacing by — \/S), and rescaling so the points again lie on a unit sphere. 
As already mentioned, the principal angles between each pair of these planes are arcsin l/y/E 
and vr/2, so d 2 c = 6/5, d 2 g = 2.6824. 

Here is another very nice packing, also found by the algorithm: for 18 planes in M 4 , use the 
binocular code consisting of the pairs (l,r), where I and r range over the vertices of a regular 
octahedron. This has d 2 = 1 (which is optimal, by the Corollary below) and d 2 = ir 2 /8 (not 
optimal). 

9 





N=2 





0,1,2 



iV=3 





N=4 





At 5 



Figure 3: Binocular codes describing best packings of N = 2, . . . , 5 planes in G(4, 2) for chordal 
distance. 

Table |2] gives the values of d 2 for the best packings we have found of iV < 50 planes in R 4 . 
Table 2: Values of d 2 for best packings found of N < 50 planes in R 4 . 



N 



6 



7 



9 



10 



11 



12 



d 2 1.5000 1.3333 1.2500 1.2000 1.1667 1.1429 1.1231 1.1111 1.0000 1.0000 



N 13 



14 



15 



16 



17 



18 



19 



20 



21 



22 



d 2 c 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9091 0.9091 0.8684 0.8629 



N 23 



24 



25 



26 



27 



28 



29 



30 



31 



32 



d 2 c 0.8451 0.8372 0.8275 0.8144 0.8056 0.8005 0.7889 0.7809 0.7760 0.7691 



N 33 



34 



35 



36 



37 



38 



39 



40 



41 



42 



d 2 c 0.7592 0.7549 0.7489 0.7477 0.7286 0.7198 0.7095 0.7066 0.6992 0.6948 



N 43 



44 



45 



46 



47 



48 



49 



50 



d? c 0.6844 0.6831 0.6809 0.6793 0.6732 0.6667 0.6667 0.6667 



Many of the other entries in Table |2j are also very beautiful: see |13[. On the other hand, 
with the exception of small values of N (6, for example), the best packings with respect to the 
geodesic distance are much uglier. That is, the packings are less symmetric, i.e. the orders of 
the automorphism groups of the best geodesic packings are almost always small, whereas those 
of the best chordal-distance packings are occasionally quite large (we will see an example in 
the next section). 
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Figure 4: Best packing of 6 planes in G(4, 2) with respect to both metrics. The left and right 
codes comprise the vertices of an icosahedron. Adjacent vertices in one code are matched with 
non-adjacent vertices in the other code. 

5. Packing n-spaces in m n 

We also applied our algorithm to look for examples of packings in E(m,n), and amassed a 
large number of examples. In looking for an analogue of Theorem 1 we collected the following 
pieces of evidence: 

(i) There is a packing of 10 planes in R 4 that forms a regular simplex. 

(ii) The 18-plane packing described in the previous section has the structure of a regular 
orthoplex (a generalized octahedron, or cross-polytope). 

(iii) We wrote a computer program to determine the lowest dimensions into which our 
library of packings in G(m, n) could be isometrically embedded. More precisely, for a given set 
of N points in G(m, n), we searched for the smallest dimension D such that there are N points 
in R D whose Euclidean distances coincide with the chordal distances between the points. 

The results were a surprise: it appeared that G(m, n) with chordal distance could be 
isometrically embedded into MP , for D = C™^ 1 ) — 1, independent of n. Furthermore the points 
representing elements of G(m, n) were observed to lie on a sphere of radius \Jn(m — n)/2m in 
R D . 

These three facts, and a number of other pieces of evidence mentioned in led us to 
the main theorem of that paper. The key idea is simply to associate to each P € G(m, n) the 
orthogonal projection map from R m to P. If A is a generator matrix for P whose rows are 
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orthogonal unit vectors, then the projection is represented by the matrix 

V = A tr A . (5) 

V is an m x m symmetric idempotent matrix, which is independent of the particular orthonor- 
mal generator matrix used to define it. Changing to a different coordinate frame in K m has 
the effect of conjugating V by an element of 0(m). With the help of @, we see that 

trace V = n . 

Thus V lies in a space of dimension ("'J 1 ) — 1. 

Let || || denote the L2-norm of a matrix: if M = (My), 1 < i,j < m, 



\M\ 



\ 



12 12 M ij = Vtrace M tr M 
i=lj=X 



For P,Q £ G(m, re), with orthonormal generator matrices A, B, and principal angles 8i, . . . , 6 n , 
an elementary calculation using (pi), (0) shows that 



2 On) 



d 2 c (P, Q) = n - (cos 2 0i H h cos 

= n - trace A tr AB tr B 

= l\\r-Q\\\ (6) 

where V, Q are the corresponding projection matrices. 

Note that if we define the "de-traced" matrix V = V — ^I m , then trace V = 0, and 
||"P|| 2 = n jT ■ We have thus established the following theorem. 

Theorem 2. The representation of n-spaces P G G(m,n) by their projection matrices V gives 
an isometric embedding ofG(m,n) into a sphere of radius \Jn{m — n)/n inR D , D = C^ 1 ) — 1, 
withd c (P,Q) = j=\\P-Q\\. 

Thus chordal distance between spaces is 1/V2 times the straight-line distance between the 
projection matrices (which explains our name for this metric). The geodesic distance between 
the spaces is l/\/2 times the geodesic distance between the projection matrices measured along 
the sphere in R D . 

Incidentally the Plucker embedding, in which members of G(m, re) are represented by points 
in projective space of dimension (™) — 1, does not give a way to realize either d c or d g as 
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Euclidean distance. Note also that the dimension of the Plucker embedding is in general much 
larger than the dimension of our embedding. 

Since we have embedded G(m,n) into a sphere of radius y/n(m — n)/m in K D , we can 
apply the Rankin bounds for spherical codes pll, and deduce: 



Corollary, (i) The simplex bound: for a packing of N n-spaces in R m , 

2 n(m - n) N 
dc ~ m 'N-l' [) 

Equality requires N < D+l = ( m ^), and occurs if and only if the N points in R D corresponding 
to the n-spaces form a regular 'equatorial' simplex, 
(ii) The orthoplex bound: for N > ( m ^" 1 )' 

, 9 n(m — n) 
m 

Equality requires N < 2D = (m — l)(m + 2), and occurs if the N points form a subset of the 
2D vertices of a regular orthoplex. If N = 2D this condition is also necessary. 

The corollary allows us to establish the optimality of hundreds of our packings. 
The case of subspaces of dimension n in E 2n or M? n+1 is especially interesting. The largest 
possible arrangements of subspaces that could achieve the two bounds are: 

m n N (simplex) N (orthoplex) 

2 1 3^ 4^ 

3 16,/ 10 

4 2 10 1/ 18 V 

5 2 15 28 ( 9 ) 

6 3 21 40 

7 3 28 ^ 54 

8 4 36 70 y/ 

Checks indicate that such a packing exists. It is known that the orthoplex bound cannot be 
achieved by 10 lines in G(3, 1), while the other cases are undecided. Our computer experiments 
strongly suggest that no set of 15 planes meets the simplex bound in G(5, 2). 

70 4-spaces in R 8 . On the other hand, it is possible to find packings of 70 points in G(8, 4) 
meeting the bound (p|). With a considerable amount of effort we determined several examples, 
of which the following is the most symmetrical. Let the coordinates be labeled oo, 0, 1, . . . , 6, 
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and take two 4-spaces generated by the vectors 

{10000000, 01000000, 00100000, 00001000} 
{11000000, 00101000, 00010001, 00000110} ' 

respectively. We obtain 70 4-spaces from these by negating any even number of coordinates, 

and/or applying the permutations (0123456), (oo0)(16)(23)(45) and (124)(365). The principal 
angles are 0,0, ^; f , j, |, \\ or ^, |, |, ^, so d% = 2, d 2 g = 7r 2 /4 (this is not even a local 
optimum with respect to geodesic distance). 

In fact, the automorphism group G of this packing acts transitively on the 70 subspaces, so 
we can obtain the packing by taking (say) the first of the spaces in ([H]) and letting the group 
act. The group has structure 2 8 ^4.g and order 5160960 (where A$ is the alternating group of 
order 8). 

6. The Miraculous Enters 

A few weeks after that packing was discovered, my colleague Peter Shor, who was studying 
fault-tolerant quantum computation, asked me about the best way to investigate a certain 
group of 8 x 8 orthogonal matrices. I replied by citing the computer algebra system MAGMA 
Hi) 0) §|) an d gave as an illustration the MAGMA commands needed to specify the group 
G of the packing of 70 4-spaces in R 8 . To our astonishment the two groups turned out (apart 
from a minor re-ordering of the coordinates) to be identical (not just isomorphic)! 

We then discovered that this group was a member of an infinite family of groups that 
played a central role in a joint paper ||] written by another colleague, A. R. Calderbank. This 
is a certain family of Clifford groups (the name is due to Wall Q, ||, f41~fl ), which may be 
constructed as follows. 

The starting point is the standard method of associating a finite orthogonal space to an 
extraspecial 2-group, as described for example in Theorem 23.10, or pq| ; Theorem 13.8. 
The end result will be the construction of various packings of re-spaces in a parent space 
V = K m , where m = 2\ As basis vectors for V we use e u , u G U = F 2 . The constructions will 
involve certain subgroups of the real orthogonal group O = 0(V, K). 

For a, b E U we define transformations X{a) £ O, Y{b) G O by 

X{a) : e u -» e u+a , Y(b) : e u -» {-lf u e u , ueU , 

where the dot indicates the usual inner product in U. Then X = (X(a) : a £ U), Y = 
(Y(b) : b G U) are elementary abelian subgroups of O of order 2 l , and E = (X,Y) C O is 
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an extraspecial 2-group 3 of order 2 2l+l (Q, Lemma 2.1). The elements of E have the form 
±X(a)Y(b), a,b &U, and satisfy 

Y(b)X(a) = (-ir b X(a)Y(b) , 

(-l) s X(a)Y(b)(-l) s 'X(a')Y(b') = (-l) a '- b+s+s ' X(a + a')Y(b + b') . 

The center of E is {±1}, and E = EfE(E) is an elementary abelian group of order 

2 2i whose elements can be denoted by X(a)Y(b), a,b £ U, where we are using the bar ~ for 
images under the homomorphism from E to E. As in Q, Theorem 23.10 we define a quadratic 
form Q : E -> F 2 by 

Q(g)-- 



if g 2 = +1 

1 ifg 2 = -I 



for g £ E, where g £ E is any preimage of g, and so Q(X(a)y(fc)) = a ■ b. 
The associated alternating bilinear form B : £7 x E — > F2 is given by 

B(gi,g2) = Q(gi + 92) + Q(gi) + Q(92) , 

for 51, 52 S and so 

S(A > (a)F(6) , A > (a / )F(6')) =a-bf + a'-b. (11) 

Then Q) is an orthogonal vector space of type r2 + (2i, 2) and maximal Witt index (cf. pq| ). 
The Clifford group L that we need is the normalizer of E in O. This has order 

2 i 2 +i+ 2 ( 2 *_i) nr^^'-i) 

(cf. Q , Section 2). For i = 3 the order is 5160960: this is the group mentioned at the end of the 
last section. L is generated by E, all permutation matrices G(A, a) 6 O : e u — ► eAu+a-, 
where A is an invertible i x i matrix over IF2 &nd a G U and. the further matrix H — (i^^ ^,), 
H u>v = 2~ i / 2 (-l) u - v , u,v£U. 

The group L acts on E by conjugation, fixing the center, and so also acts on E. In fact L 
acts on E as the orthogonal group + (2i,2) (||, Lemma 2.14). 

This Clifford group L has arisen in several different contexts, providing a link between the 

the problem of packing in Grassmannian spaces, the Barnes- Wall lattices (see 0], @, |p^| , 

the construction of orthogonal spreads and Kerdock sets ||] , and the construction of quantum 
3 |28|, p. 349 
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error-correcting codes |J, |fL2]1 . It also occurs in several purely group-theoretic contexts - see 
H for references. 

The connection with quantum computing arises because if certain conditions are satisfied 
the invariant subspaces mentioned in Theorem 1 form good quantum-error-correcting codes 

0, @. 

7. The construction from totally singular subspaces 



We give one theorem, taken from |10fl , to illustrate how the Clifford group may be used to 
construct packings in Grassmannian spaces. 

A subspace S C E is totally singular if Q(g) = for all g G S. Then dim S < i, and if 



dim S = i then S is maximally totally singular. It follows from (11) that the preimage T C E 
of a maximally totally singular space T is an abelian subgroup of E, of order 2 4+1 . T contains 
—J, and has 2 4+1 linear characters, associated with 2 l mutually perpendicular 1-dimensional 
invariant subspaces forming a coordinate frame J-(T) C V (||, Lemma 3.3). 

Since L acts as + (2i,2) on JE, L takes any ordered pair of maximally totally singular 
subspaces that meet in {0} to X and Y respectively. The corresponding coordinate frames in 
V are 

HX) = {el = ^-Y. (-^^u :veU} (12) 

and 

F{Y) = {e u : u G U} , (13) 

respectively. 

If S C T has dimension fe, its preimage S C £" has 2 fc+1 linear characters, and 2 fc distinct 
invariant subspaces, each of which is spanned by 2 % ~ k of the vectors in J-(T). 

The following theorem produces many good Grassmannian packings. For the proof see [ 10 1 . 



Theorem 3. Given k, with < k < % — 1, the set of all invariant subspaces of the preimages S 
of all (i — k)- dimensional totally singular subspaces S of E is a packing of N planes in G(2 l , 2 k ) 
with minimal distance d = 2^ k ~ 1 '' 2 ) where 

i-l 



N 



■\i—k 



k 



ri( 2j + i ) 

j=k 



16 



and 

' i ] _ (2* - 1) . . . (2 i ~ k+1 - 1) 
k \ ~ (2 fe -l)...(2-l) 

is a Gaussian binomial coefficient. 

Examples. Taking k = in the theorem we obtain a packing of 

(2 + 2)(2 2 + 2)...(2* + 2) 



lines in G(2 l , 1) with minimal angle 7r/4 (as in [|37|). These are the lines defined by the minimal 
vectors in the 2 l -dimensional Barnes- Wall lattice together with their images under H (cf. [ji~4"f , 
p. 151). 

With i = 2, k = 1 and i = 3, fc = 2 we obtain two important special cases: 18 points in 
Cr(4, 2) and 70 points in G(8,4). More generally, when k = i — 1 we obtain the packing of 

/(i) = 2(2 i -l)(2^ 1 + l) 



points in G(2 i ,2 i ~ 1 ) with d 2 = 2 l ~ 2 that is the main result of j37|. These packings meet the 
orthoplex bound of (||) and are therefore optimal. An explicit recursive construction for the 
special case k = i — 1 is given in 



For k = 1 and k = i — 2 we obtain two further sequences of packings whose existence was 
conjectured in jj7|]. 

The construction given in the theorem can be restated in an equivalent but more explicit 
way as follows. Let Po be the 2 fc -dimensional space spanned by the coordinate vectors e u , 
where u G U is of the form 00 ... * ... *, with i — k initial zeros. Then the packing consists 
of all the images of Pq under the group L. 



The paper |1C] also contains many other Grassmannian packings. 

Space does not permit any discussion of the application of the Clifford group to construct 
quantum error correcting codes: for this see [[TJ]. 

We conclude by mentioning the present status of the existence of the packings listed in 
Eq. (^). No packing of 10 points in G(3, 1) meeting the orthoplex bound can exist, and we 
conjecture no packing of 15 points in G(5, 2) meeting the simplex bound can exist. We have 
found an infinite family of packings of p(p + l)/2 points in G(p, (p— l)/2) meeting the simplex 
bound, generalizing the packing of 28 points in G(7, 3) |l(J. As mentioned above, we have also 
found an infinite family of packings of m 2 + m — 2 points in G(2 m , 2 m ~ 1 ) meeting the orthoplex 
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bound, generalizing the packings of 18 points in G(4, 2) and 80 points in G(8,4). We would 
very much like to know if packings of 21 or 40 points in G(6, 3), 54 points in G(7, 3) or 36 
points in G(8, 4) meeting the bounds exist. 

Acknowledgements. I should like to express my thanks to my coauthors (Rob Calderbank, 
John Conway, Ron Hardin, Eric Rains and Peter Shor) of the papers to which this article 
serves as introduction. 
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